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We report on strong renormalization encountered in periodically driven interacting quantum dots 
in the non-adiabatic regime. Correlations between lead and dot electrons enhance or suppress 
the amplitude of driving depending on the sign of the interaction. Employing a newly developed 
flexible renormalization group based approach for periodic driving to an interacting resonant level 
we show analytically that the magnitude of this effect follows a power law. Our setup can act as a 
non-Markovian, single-parameter quantum pump. 
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Lately, signal and information processing at the quan¬ 
tum level has attracted sizable interest. It can be realized 
in quantum dot geometries which can be manufactured 
very precisely, e.g., in semi-conductor heterostructures. 
These provide the basic building blocks of future nano- 
electronic and quantum information devices [T]. An ex¬ 
perimentally routinely implemented non-equilibrium dot 
setup is the one of periodically varying external fields 
[a [3]. utilizing as an additional parameter the frequency 
Ul of the field one can achieve novel control of the quan¬ 
tum dot’s transport properties. Examples of such include 
a vanishing of certain tunneling amplitudes of the dot 
structure, which is analogous to the coherent destruction 
in two-level systems mis], and pumping of an integer 
number n of the elementary charge e during one cycle 
EHH]. The latter is important for metrology, where the 
quantized charge ne could replace the currently used elec¬ 
tric current standards [5]. 

In the adiabatic limit of a small rate of change of the 
external field the latter can be used as a small parameter 
in analytic calculations. This was employed to gain valu¬ 
able insights into the physics of periodically driven quan¬ 
tum dots [IMS]. The regime in which Vl is of the order 
of the (average) tunneling rate T or larger, however, was 
treated less extensively [TMTH] . but showed interesting 
physics, e.g. one parameter pumping [l9l 1^ . In partic¬ 
ular, describing the many-electron interplay encountered 
in non-adiabatically driven quantum dots in presence of 
correlations remains a formidable challenge iniiiijiig. 
We here report on progress on this. 

Following the usual procedure [lUHiHi uni na , we 
model the influence of the external time-periodic fields by 
time-periodic parameters p{t) of the quantum dot setup. 
In general, all Hamiltonian parameters, such as the dots 
level positions and the tunneling amplitudes might ac¬ 
quire a time dependency. Given these one can compute 
the relevant observables, such as the current through the 
quantum dot, e.g. using scattering theory [TUI [TH] - We 
emphasize that depending on the setup under consid¬ 
eration obtaining closed analytical expressions for ob¬ 
servables of this time-dependent single-particle problem 



FIG. 1. (Color online) A periodically varied tunneling ampli¬ 
tude between the left reservoir and the quantum dot Ti,{t) is 
strongly renormalized by an interaction U. The upper panel 
shows tl (t) for one period, which is designed to yield a renor¬ 
malized hopping amplitude of rectangular line shape (see the 
lower panel). It turns out to dependent on sgn(?7). The cen¬ 
tral panel depicts the studied model, which consists of a quan¬ 
tum dot coupled to two electronic reservoirs via hopping am¬ 
plitude r and a dimensionless interaction U. For U = TO.S/tt 
the emergent universal energy scale is Tk{Uti — —0.5)/D = 
4.41 • 10~'^ and Tk{U-k = 0.5)/D = 2.865 • 10'^. The Fourier 
components of the hopping are IrL^fc/rol = QA'Kk!'~^^^\ As 
shown in the lower panel the driving is rectified for positive 
U (red lines), while it is amplified for negative U (blue lines), 
where both effects enhance with increasing driving frequency 
D following a power law. 


might turn out to be rather involved (for examples rel¬ 
evant here see below and [53], [55]). In the presence of 
two-particle interactions a many-particle method must 
be employed. We use a newly developed extension of the 
functional renormalization group (ERG) [26] to periodi¬ 
cally driven systems applicable at arbitrary driving fre¬ 
quency. We show that the emergent correlations lead to 
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a strong, qualitative renormalization of the time-periodic 
parameters p{t) —)■ in the anti-adiabatic regime 

n ^ r. This affects not only their amplitude but their 
entire line shape (see Fig. [^. For typical quantum dot 
models the renormalization is not accessible by pertur¬ 
bation theory in either the two-particle interaction or the 
tunnel coupling. It can be exploited to flexibly design the 
effective driving and thus the time periodic current by 
controlling the frequency (see Fig. [^. Within our trun¬ 
cation of the hierarchy of FRG flow equations the p"'®"(t) 
can directly be plugged into the non-interacting expres¬ 
sions for observables. Following this two step procedure 
of (i) setting up and solving RG equations for the time- 
dependent parameters and (ii) substituting them into 
noninteracting expressions we device a method which al¬ 
lows one to obtain numerical results for arbitrary dot 
parameters and analytical ones in limiting cases such as 
e.g. the anti-adiabatic regime of interest here. 

To depict the parameter renormalization in a trans¬ 
parent fashion, we concentrate on a spin polarized, single 
level quantum dot. It hosts a fermionic degree of free¬ 
dom and is coupled to two infinite, equilibrated reser¬ 
voirs. We consider the left hopping Ti,{t) = J2k 
to be time-periodic, allowing for an arbitrary line shape, 
i.e. arbitrary Fourier amplitudes tl.^. The right hopping 
is set constant TR(t) = = R); where we focus on left- 

right symmetric mean values. Additionally, we consider 
a capacitive density-density type interaction U between 
electrons residing in the reservoirs and on the dot as well 
as an onsite energy e, which for simplicity are both as¬ 
sumed to be time-independent. The underlying model 
is known as the interacting resonant level model (IRLM) 
[?7l - l55| , which is the standard model for correlated quan¬ 
tum dots dominated by charge fluctuations. The Hamil¬ 
tonian is H = Hdot + Hcoup + Hres, with iJdot = ed'^d 
and 


-^coup — ^ ^ ^a(^) 

q',Q:=L,R 


+ idU-- 




H : • ( 1 ) 

q,q^ ,a—L,R 


We use second quantization notation, where the dot and 
reservoir operators are d and c, respectively, and : ... : 
denotes normal ordering. The reservoirs given by iilres 
only enter via their local density of states Pa{(^)- For 
convenience we choose a symmetric Lorentzian Pai<-o) = 
Pa^ with pa^ = ^, where the bandwidth D ^ 

|rQ,(t)|, |e|. In this (scaling) limit the details of the energy 
dependence of pa are irrelevant. We assume Ua = U but 
consider both cases of positive and negative interaction. 
The latter might effectively be realized for a quantum 
dot coupled to phonons, with the phonon frequency in 
the adiabatic limit [M] . 

The IRLM can be mapped to the anisotropic Kondo 
model [35] and our results are of relevance for this model 


as well. The Kondo effect under periodic driving was 
studied earlier [22l [36l [37] employing the single-impurity 
Anderson and the isotropic Kondo model. 

We show that in anti-adiabatic driving the correla¬ 
tions induced even by a small U, strongly renormalize 
the Fourier coefficients ru.fc —>■ This is captured by 

the central analytical result of our Letter |21| 


--f{T0,U)T0, (2) 

~ L!-^(*fc)-V,fe, k^O (3) 

where H is the largest scale (besides the bandwidth D). 
We furthermore assume a small amplitude of driving 
[tl j.^o/'^’ol ^ 1- l^ote, that the formalism we device can 
also be used beyond these limits and for more general se¬ 
tups with all dot parameters being time-periodic. While 
'^L™=o @ is independent of Q, the higher harmon¬ 

ics Eq. ^ are suppressed or amplified in a power-law 
fashion with increasing driving frequency, depend¬ 
ing on the sign of the interaction U. The factor k~^ 
leads to a change of the line shape of T£®"(t) after trans¬ 
forming back to time space. Gonsequently, varying H for 
a fixed C/, one finds a constant mean value and equal 
line shape of the renormalized r£®"(t), but a strongly in¬ 
creased or suppressed amplitude. This renormalization 
can be exploited as an amplifier or rectifier. Let us as¬ 
sume that one is interested in achieving a given line shape 
of the hopping amplitude and thus a given shape 

of the time-periodic current. We here choose the former 
to be rectangular, but any form could be targeted. At 
a given value of U we design the external driving such 
that {ik)~^Tj^ ^ is equal to the Fourier coefficients of a 
periodic rectangular pulse. This is depicted in Fig. for 
[7 ^ 0 at e = 0. Here we already use the emergent scale 
Tk as the unit of energy, which will be introduced next. 
Power laws in the IRLM —Already in equilibrium observ¬ 
ables of the IRLM, e.g. the charge susceptibility 



e=0 




(4) 


exhibit power-law scaling with interaction dependent ex¬ 
ponents. Plain perturbation theory (in U) thus leads to 
logarithmic terms. Various approaches were devised to 
resum these The susceptibility can be used to 

define an emergent low-energy scale [28], Tk = —2/ (-Try). 
Universal power laws were encountered in the (static) 
bias voltage driven non-equilibrium steady state, e.g. in 
the current-voltage characteristics US] [5TH551155] . Fur¬ 
thermore, even the transient relaxation dynamics shows 
power laws in the time variable with interaction depen¬ 
dent exponent [33l |39l EO]- This renders the dynamics 
of the periodically driven IRLM a likely candidate for 
power-law renormalization. 

Renormalization group approach to periodically driven 
quantum dots —Here we introduce a reformulation of the 
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general time-dependent FRG of Ref. [33] to the time- 
periodic problem of interest by making use of the Flo- 
quet theorem. This allows one to efficiently treat time- 
periodic setups with arbitrary and gain analytical in¬ 
sights in limiting cases, which are difficult to obtain oth¬ 
erwise. Within the FRG framework a suitable chosen 
cutoff A is introduced. The quantum many-body prob¬ 
lem is rephrased in terms of a coupled set of infinitely 
many differential flow equations with respect to A for 
the one-particle irreducible (Keldysh) vertex functions 
[26| . The only approximation consists in truncating this 
set, here by suppressing the renormalization of the two- 
particle vertex (effective interaction), which is 0(C/^). As 
a consequence, only the self-energy acquires a RG flow. 
The effect of the interaction U is incorporated succes¬ 
sively from high to low energies, resulting in renormalized 
single-particle parameters (self-energy). 

We are interested in the periodically driven steady 
state and thus assume that the driving as well as the 
dot-reservoir-coupling were switched on in the distant 
past. As usual in non-equilibrium retarded, advanced 
and Keldysh Green’s functions are introduced and writ¬ 
ten in standard matrix form m 

( ^ret (^K \ 

0 (Jadv j ■ (5) 


An analogous structure is found for the self-energy. We 
employ the auxiliary infinite-temperature reservoir cutoff 
described in [42] . to introduce the flow parameter A in 
every non-interacting Green’s function. The two-times 
Green’s functions are determined by the Dyson-equation 



eeff(t) 




( 6 ) 


where the single-particle part of the Hamiltonian eeff(t) 
incorporates the effect of the self-energy; the latter re¬ 
mains time local in our approximation |42j . This dif¬ 
ferential equation can be solved for time-periodic eeff(t) 
exploiting the Floquet theorem [13]. One writes 

G^{t,uj)= (7) 


and 

Gf(cc) = g/dte*'=^*G^(t,a;), (8) 

with X G {ret,K,adv}, fc G Z being a discrete Flo¬ 
quet index, and w a continuous real frequency |33fl46j . 
The self-energy and observables of interest are trans¬ 
formed accordingly. Single-particle-like indices are easily 
incorporated in a FRG framework increasing the com¬ 
plexity only in a polynomial fashion |47j . Formulated 
in this way the FRG approach to the steady state of 
the time-periodic problem is closely related to the static 
non-equilibrium FRG [33]. The technical details of the 


Floquet FRG including the full RG flow equations for 
arbitrary (amplitude and line shape) time-periodic dot 
parameters will be presented elsewhere |48| . For gen¬ 
eral parameters the flow equations must be solved nu¬ 
merically with the only additional challenge compared to 
established equilibrium and static non-equilibrium FRG 
schemes that a sufficient number of higher harmonics k 
must be kept; this can be chosen to be G(IOOO) without 
exceeding available numerical resources. We note that all 
our figures show the results of the numerical solution of 
the full set of flow equations. 

In the limit of small amplitudes Ta^k one can derive 
simplified RG equations which for large can then be 
solved analytically. This allows us to gain valuable an¬ 
alytical insights into the mechanism of power-law sup¬ 
pression/amplification described above. For simplicity 
we concentrate on e = 0. In this case the flow equations 
of the hopping amplitudes are the only relevant ones, 
while the renormalization of the onsite energy is zero 
(or G([/^) for e ^ 0) [33]. To leading order in and 
U the differential flow equations for k decouple and 
with 'rl,k=o = ^R.fc=o = '^0^ are 

^ A ^ _rr _ 

ATL,k^0 A^ I /.f Hop I ikQ\A , 2i\To\'^kQ , 4|toP 

D > D D* 

(9) 

and 

U— _ . (10) 

{A/D)^+A/D + 2{to/D)^ ^ ’ 

As emphasized above the flow of Tq is independent of ft 
while in the one of k^o^ acts as a cutoff. One can 
solve these equations in a closed form (see [33]) yielding 
Eqs. ^ and ([^ for large H. The consequences of this 
renormalization as well as its potential as an amplifier or 
rectifier, were discussed above. 

Observables —In our truncation of the exact hierarchy of 
flow equations one can substitute the renormalized time- 
periodic parameters into the non-interacting expressions 
of observables, e.g. the current through the quantum 
dot structure or the occupancy of the dot. This allows 
one to obtain (approximate) closed form expressions in 
the presence of interactions and is a particular advantage 
of our approximate RG formalism. We here focus on 
the time-periodic current Ja{t) = leaving 

reservoir a which is given by 






— G'fl_k/{uj + +[—k^k] 


( 11 ) 


where the mean value is the zeroth component defined 
as Ja,k=o = ^ E'7a(i) and '£^ak denotes the self-energy 
of the reservoir a [33]. The components of the Green’s 
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FIG. 2. (Color online) The linear conductance Ghit) for 
the non-interacting case (upper panel) as well as renormal¬ 
ized by the interaction (lower panel) for two values of driv¬ 
ing frequency for the parameters of Figj^ The amplifica¬ 
tion/rectification effect of the hopping is reflected. The inset 
shows that the (odd) higher harmonics (even ones are zero) 
of the conductance Gl,*: collapse for different k if Gl,*: and 
the driving frequency are scaled by a factor of and k, 

respectively. 



FIG. 3. (Color online) The upper panel shows the pumping 
power JL,fc=o(f2) in units of JL,fc=o(f^ = Tk) for a single pa¬ 
rameter pump, tl = To -I- Ar sin(flt) with and ^ = 0.05. 
It shows again an amplification for negative as wel? as a rec¬ 
tification for positive values of U and inherits the power-law 
behavior of the hopping amplitude. The lower panel shows 
the logarithmic derivative of JL,fc=o indicating the effective ex¬ 
ponent a of the power law behavior, which is in a very good 
agreement with the predicted exponent (dotted lines). The 
Tk/D are {4.4M0“'^, 8.4•10■^ 2.2-10“®, 5.55-10"®, 1.24-10"®, 
2.5 • 10“®, 4.62 • 10"®, 8.43 • 10"®, 1.28 • 10"^, 1.96 • 10"^, 2.865 • 
10"^} for Utv = {-0.5, -0.4,..., 0.,..., 0.4, 0.5}. 


function and self-energy contain the interaction via the 
renormalized parameters. We exemplify the impact of 
the parameter renormalization of the periodically driven 
dot on the current for two examples. 

Linear response current under the driving of Fig. [^-We 
first consider the driving protocol of Fig. and study the 
linear response of J{t) to an additional small static bias 
voltage V = Hl — Mr between the left and right reservoir, 
where V is much smaller than any other scale. The effec¬ 
tive single-particle problem of a dot with a time-periodic 
left hopping of rectangular shape (see Fig. is rather 
involved as the different Floquet ‘channels’ indicated by 
k are all coupled. We thus refrain from further simpli¬ 
fying Eq. ( [TT] ) . We define Gq (t) = limv_>o dJa (t) /dV as 
the linear conductance and analogously its Fourier coef¬ 
ficients Gak = hmy_>o dda fc/dE. The upper panel of 
Fig I shows the bare, i.e. non-interacting Gi^{t) resulting 
from the external driving of Fig. for LI/Tk = 10 and 
100. As seen in the lower panel the overall amplitude of 
the renormalized conductance is also rectified/amplified 
as expected from Eq. 

The inset of Fig. [^depicts the (odd) higher harmonics 
of the conductance (even ones are zero) for different U. 
A collapse of the conductance for different k is found at 
fixed U when rescaling Gl.A: by k^~^ and 17 by /c. This 
shows that although Gl( 7) depends on the renormalized 
dot parameters in an involved way, universal scaling es¬ 
tablished in the IRLM in and out-of-equilibrium is also 
realized for periodic driving. 

Pumping —Finally, we consider the quantum dot as a 
single-parameter pump. In quantum pumping, one usu¬ 
ally studies a sinusoidal time-dependence of the dot pa¬ 


rameter(s) [TUHIS]- Assuming Ti,(t) = tq + ATsin(I77) 
we follow this convention. By breaking the particle- 
hole symmetry via (arbitrarily) choosing e = O.dTK, 
J'L,k=o 7 ^ 0 for non-adiabatic driving. JL,fc=o describes 
the average amount of pumped charge per unit time and 
thus indicates the pumping power. It can be related to 
the pumped charge per cycle Q = rTL,fe=o [iQl [HI HI]- 
Following our two step approach we derived [24] a closed 
form expression for JL,fc=o based on Eq. (Ill (see also 
[25] ) - which is straightforward as the Floquet channels 
are only ‘weakly’ coupled for sinusoidal driving - and 
substitute the renormalized parameters, which yields to 
order (|T£®f |/t£®q )^ and for large 17 


J n»e 

UL,fc=0 — 



ren 

L,0 


2 

Tk arctan 



( 12 ) 


In the upper panel of Fig. |^we show JL,fe=o for different 
[7. A clear power law Jh,k=o ~ with 7 = —2U man¬ 
ifests in the pumped charge at large driving frequency, 
which is a direct consequence of the renormalization of 
the tunneling amplitude r£°" Eq. (|^. Our main result of 
non-perturbative power-law amplification or rectification 
can thus be observed in the straightforwardly measur¬ 
able pumping power. We note that a numerical solution 
of the full set of FRG flow equations is sufficiently accu¬ 
rate to extract the exponent 7 via a logarithmic deriva¬ 
tive dln(JL,fc=o)/'^ln(^^)j (implemented as centered dif¬ 
ferences) which is a numerically very sensitive measure; 
see the lower panel of Fig. This exemplihes that from 
the data for observables detailed information can be ex- 
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tracted which is crucial in cases in which the effective 
scattering problem is to involved to gain analytical in¬ 
sights. 

The pumped current is a consequence of the non- 
Markovian nature of the reservoirs and it is not obvious 
how to capture it in standard approaches relying on a 
tunneling rate which depends on a single time argument. 
In contrast to the quantum pump described in Ref. ^D] 
(realized for ^ T) it cannot be understood in terms of 
interaction induced effective phase differences of the dot 
parameters. 

Conclusion —We have shown that in quantum dots, 
where one hopping amplitude between the quantum dot 
and the electronic reservoir is varied periodically, the ef¬ 
fective hopping amplitude is renormalized strongly by 
correlations, resulting in power-law scaling of the Fourier 
coefficients of the hopping. This is reflected in the line 
shape and the overall amplitude of the renormalized hop¬ 
ping. Depending on the sign of the interaction a power- 
law amplification or rectification is achieved. This renor¬ 
malization in turn manifests in observables such as the 
linear conductance in the presence of an additional in¬ 
finitesimal static bias voltage or the pumped charge of a 
single parameter quantum pump. The particular setup of 
only one hopping amplitude chosen as time-periodic can 
be easily generalized to one with any one (as e.g. e(t)) 
or even more parameters being time-periodic within our 
newly developed RG approach. The renormalization ef¬ 
fects are relevant to experiments, as they dominate for 
driving frequencies of the order of IOTk, which for dot 
geometries with an effective hybridization of 1 MHz is 
within an accessible regime. 
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Renormalization in periodically driven qnantnm dots 
Supplementary Material 


Here we provide (a) the analytical solution of the flow Eqs. (9) and (10) as well as (b) the derivation of the 
closed expression of the pumping power of the effective non-interacting model given in Eq. (12) of the main text. 
This exemplifies how we are able to gain analytic insights into the renormalization effects employing our functional 
renormalization group scheme. After the flow equations have been derived these can be solved analytically (in certain 
limits). Due to the truncation used here it suffices to plug the renormalized parameters obtained from the flow 
equations into the non-interacting expressions of the observables of interest. 


(a) Analytical solution of the flow equations 

The differential Eqs. (9) and (10) can be solved analytically, yielding = tq) 
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In the last step, we have specified that the frequency fl which yields Eq. (3) of the main text (disregarding the 

irrelevant proportionality ~ D^^). 


(b) Single parameter quantum pump 

It is a particular advantage of our functional renormalization group scheme that once the flow equations have 
been solved observables can be calculated straightforwardly by considering an effectively non-interacting system with 
renormalized parameters. Those renormalized parameters include (via the RG procedure) the influence of the two- 
particle interaction. This allows for a derivation of analytic expressions for the observables, which depend on the 
details of the driving protocol. Here we consider as an example Ti^{t) = tq -(- Arsin(nt) and finite e, which includes 
the quantum pump described in the main text. We aim at a closed form expression for the pumping power Jl,/:=o 
to order C9(r£ j^) in the high frequency regime. Given the renormalization of tq and (see Eqs. (SI) and (S2) or 
Eqs. (2) and (3) of the main text) we this way provide analytical insights into the power-law renormalization of Jl,/:=o 
discussed in the main text. 

It is straightforward to solve the scattering problem of non-interacting fermions for the resonant level model with 
renormalized r^.k numerically. However, to reveal insights into the pumping mechanism (and its power-law suppres¬ 
sion/amplification) we choose a more instructive path in the following. Remind that the reservoirs are modeled via 
a Lorentzian density of states pa{oj) = with D being the largest scale (scaling limit). This description is 



2 


k 

A 


-J T. ^ 




c 

c 

c 


FIG. Si: (Color online) The non-interacting driven quantum dot can be formally mapped to a non-driven one-dimensional 
system by interpreting the Floquet index fc as a spatial one. 
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FIG. S2: (Color online) To 0 (t£ i) the pumped current can be evaluated in an effective three-terminal setup. To the left and 
right the energy dependent density of states of the three reservoirs are shown. 


equivalent to keeping one site of the reservoir explicit and coupling it to a structureless reservoir with hybridization 
D. In this language the single reservoir site kept explicitly is coupled to the central level via tq. Next we note that 
including the Floquet index k for the periodically driven system is mathematically identical to enlarging the dimen¬ 
sionality of the system (the index k then takes the role of an additional spatial index, see [25]). Along the infinite 
‘fc-direction’ replicas of the initial system with all energy levels shifted by kQ, have to be included. The different 
replicas are coupled via TL,k (where the Floquet index k indicates the range of the coupling in the ‘fc-direction’). 
With this construction one can map the non-interacting periodically driven zero-dimensional quantum dot to an one¬ 
dimensional system in the steady state (no periodic driving). As we consider a sinusoidal form of only tl,i is 
non-zero and couples neighboring replicas. The mapping is sketched in Fig. SI 

We aim at an analytical expression for the current J through the central channel A: = 0 of the mapped model 
shown in Fig. SI to C>(r| j^) in the large frequency fl limit. To C>(r| j^) the influence of the replicas below k = —1 and 
above k = 1 (dotted in Fig. SI) can be neglected and the contribution of the fc 7^ 0 channels to the k = Q one can be 
evaluated in their respective equilibrium (where all i are set to zero). 

To sum up, the scattering problem through the central channel A: = 0 is thus dressed by temporary excursions of 
particles into the A: = ±1 channel (and back). The temporary excursions exhibit an equilibrium background in the 
k = ±1 channel. We aim at the current leaving the left reservoir [denoted by J(= JL,fc=o) in Fig- SI]. Left and right 
reservoirs are held at the same chemical potential and thus to (!I('r£ j^) there is no net current JL,fc=o- Also temporary 
excursions of particles from the central site at A: = 0 to the left site of the A: = ±1 channel and back to the central 
site of A: = 0 always respect left-right symmetry and give no net contribution to the current Jh,k=o from left to right. 
Particle excursions from the central site at A: = 0 to the left site of the A: = ±1 channel and back to left site of A: = 0 
break the left-right symmetry, but are suppressed by 1/Q, and provide no contribution for large frequency. The only 
relevant process to j^) is thus the one from the left site of A: = 0 to the central site of A: = ±1 and back. Now we 
consider the left site of the k = 0 channel, which within A: = 0 is coupled much more strongly to the left (~ D) than 
to the right (~ tg/D). Therefore the infiuence of the coupling of the left site to the central site of the k = 0 channel 
will be neglected in leading 0{1/D). Calculating J-L^k=o thus effectively maps to the problem of a single site (former 
left site of A: = 0) coupled to three reservoirs. This is depicted in Fig. S2. 
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Using the Landau-Biittiker equation [49] one now finds Eq. (12) of the main text by 
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for e ^ n <C -D and where Tk = Tr^jD. An expansion in the limit of small driving frequency 52 yields the leading 
contribution ~ e 522 renders an adiabatic approach insufficient to describe the pumped charge and thus underlines 
the relevance of the non-adiabatic physics at play. In the opposite limit of large driving frequency 52 2> e, plugging 
in the renormalized values for tq and t^ i from Eqs. (SI) and (S2) where only the second depends on 52 one recovers 
the power-law JL,fc=o ~ 52 ^ 2 C^ mentioned in the main text (compare also Fig. 3 of the main text). Note that Eq. (S3) 
has the structure of a finite bias current through a non-interacting resonant level with bias voltage replaced by the 
onsite energy t. 



